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We report on a study of three-body charmless decays B'^ K^h'^h~ based on a 29.1 fb~' 
data sample collected with the Belle detector. With no assumptions on the intermediate mech- 
anisms, the following three-body branching fractions have been measured for the first time: 
B{B+ K+-K--K+) = (55.6±5.8±7.7) x 10"'* and B{B+ ^ K+K'K+) = (35.3±3.7±4.5) x 10"^ 
We present the first observation of the decay B'^ /o(980)ii'^ with a branching fraction prod- 
uct of B{B+ fo{mO)K+) X S(/o(980) ^ Tr+Tr") = {9.6tlit\itli) x lO"''. This is the first 
reported example of a B meson decay to a scalar pseudoscalar final state. We also report the first 
observation of B+ K*{892)°tt+ decay with a branching fraction of B{B+ K*{892)°tt+) = 



PACS numbers: 13.20.He, 13.25.Hw, 13.30.Eg, 14.40.Nd 



I. INTRODUCTION 



During the last few years, a considerable amount of new information on charmless hadronic decays of B mesons has 
been reported, primarily by the CLEO Collaboration. The discoveries of the B Kn and B —>■ nn decay modes 0| 
have provided a real basis for searches for direct CP- violating effects in the B meson system. 

However, because of large combinatoric backgrounds, studies of charmless B decays have concentrated mainly on 
two-body decay processes. Three-body decays could significantly broaden the study of B meson decay mechanisms 
and provide additional possibilities for direct CP violation searches. In this paper, we report the results of a study 
of charged B meson decays to three charged particle final states Kmr, KKt:, and KKK, where no assumptions are 
made about intermediate hadronic resonances. We also present the results of a study of quasi- two- body intermediate 
states in the if+Tr+Tr" and K^K'^K~ final states. The inclusion of charge conjugate states is implicit throughout 
this work. 

The data sample used for this analysis was collected with the Belle detector operating at the KEKB asymmetric 
energy e+e" collider It consists of 29.1 fb"' taken at the T(4S') resonance, corresponding to 31.3 x 10^ produced 
BB pairs, and 2.3 fb~' taken 40 MeV below the BB production threshold to perform systematic studies of the 
e+e" qq background. 



II. THE BELLE DETECTOR 

The Belle detector js) is a large-solid- angle spectrometer based on a 1.5 T superconducting solenoid magnet. Charged 
particle tracking is provided by a three layer double-sided silicon vertex detector (SVD) and a 50 layer cylindrical drift 
chamber (CDC) that surround the interaction region. The charged particle acceptance covers the laboratory polar 
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angle between = 17° and 150°, corresponding to about 92% of the full solid angle in the center-of-mass (cm.) frame. 
The momentum resolution is determined from cosmic rays and e+e" — > n~ events to be (Jpjpt = (0.300 0.19pt)%, 
where pt is the transverse momentum in GeV/c. 

Charged hadron identification is provided by dE/dx measurements in the CDC, an array of 1188 aerogel Cerenkov 
counters (ACC), and a barrel- like array of 128 time-of-flight scintillation counters (TOF). At large momenta 
(> 2.5 GeV/c) only the ACC and dE/dx are used for separation of charged pions and kaons since here the TOF 
provides no additional discrimination. 

Electromagnetic showering particles are detected in an array of 8736 CsI(Tl) crystals that is located in the magnetic 
volume and covers the same solid angle as the charged particle tracking system. The energy resolution for electromag- 
netic showers is (Je/E = (1.3® 0.07/i?ffi0.8/i?^/'')%, where E is in GeV. Electron identification in Belle is based on a 
combination of dE/dx measurements in the CDC, the response of the ACC, and the position, shape and total energy 
deposition (i.e. E/p) of the shower registered in the calorimeter. The electron identification efficiency is greater than 
92% for tracks with pi^h > 1-0 GeV/c and the hadron misidentification probability is below 0.3%. 

The magnetic field is returned via an iron yoke that is instrumented to detect muons and Kl mesons. We use a 
Monte Carlo simulation to model the response of the detector and determine acceptance Q. 

III. EVENT SELECTION 

Charged tracks are required to satisfy a set of track quality requirements based on the average hit residual and 
on the distances of closest approach to the interaction point in the plane perpendicular to the beam and the plane 
containing the beam and the track. We also require that the transverse track momenta be greater than 0.1 GeV/c to 
reduce the low momentum combinatoric background. 

Charged kaon candidate tracks are selected with a set of FID criteria that has about 90% efficiency, a charged 
pion misidentification probability of about 8%, and a negligible contamination from protons. We also reject tracks 
that are identified as electrons. Since the muon identification efficiency and fake rate vary significantly with the track 
momentum, we do not reject muons to avoid additional systematic error. 

We reconstruct B mesons in three charged track final states with at least one positively identified kaon. The 
candidate events are identified by their cm. energy difference, AE = (^^ Ei) — E^,, and the beam constrained mass, 
Mbc = \/E'^ — iJ2i'Pi)^j where E^ = •\/s/2 is the beam energy in the cm. frame, and pi and Ei are the cm. three- 
momenta and energies of the candidate B meson decay products. We select events with Mbc > 5.20 GeV/c^ and 
|A£:| < 0.20 GeV, and define a signal region of |Mbc - Mb\ < 9 MeV/c^ and |A£:| < 0.04 GeV and two AE sideband 
regions defined as -0.08 GeV < A£; < -0.05 GeV and 0.05 GeV < AE < 0.15 GeV |||. The selection of sideband 
regions was based on a Monte Carlo study and was done in such a way that the relative fraction of BB and qq events 
match that of the signal region. For the normalization factor between the sideband and signal data samples we use 
the ratio of areas, namely 0.62. 

To evaluate signal and background levels, we require that one of AE or Mbc fall in its signal region and examine 
the distribution of candidates in the other, fitting to the sum of a signal distribution and an empirical background. 
The signal shapes in Mbc and AE distributions are parameterized by a Gaussian and sum of two Gaussians with 
the same mean, respectively. The width of the Mbc distribution is primarily due to the cm. energy spread and is 
expected to be the same for each channel; in the fit we fix it at the value UA/b^ = 2.9 MeV/c^ determined from 
the B'^ ZJ'^TT+^DO K^TT^ events in the same data sample. The AE shape for the signal is also determined 
from the B~^ — > D^n'^ events. For the Mbc projection, we parameterize the background with the empirical function 
/(Mbc) ex. y/l — x^ exp[— ^(1 — x^)], where x = Mhc/Eh and ^ is a parameter Q. We fix the ^ value from a study of 
data below the BB production threshold. We represent the AE background shape with a linear function and restrict 
the fit to the range -0.1 GeV < AE < 0.2 GeV g. 

IV. BACKGROUND SUPPRESSION 

An important issue for this analysis is the suppression of the large combinatoric background which is dominated 
by e+e~ qq continuum events. We suppress this background with variables that characterize the event topology. 

Since the two B mesons produced from T{AS) decay are nearly at rest in the cm. frame, the angles of the decay 
products of the two i?'s are uncorrelated and the events tend to be spherical. In contrast, hadrons from continuum qq 
events tend to exhibit a two-jet structure. We use Othi, which is the angle between the thrust axis of the B candidate 
and that of the rest of the event to discriminate between the two cases. The distribution of | cos^thrl is strongly peaked 
near |cos0thr| = 1-0 for qq events and is nearly flat for BB events. We require jcos^thrl < 0.80 for all three-body 
final states; this eliminates 83% of the continuum background and retains 79% of the signal events. 
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After imposing the cos^thr requirement, the remaining qq and BB events still have some differences in topology that 
are exploited for further continuum suppression. We divide the space around the B candidate thrust axis into nine 
polar angle intervals of 10° each; the i-th interval covers angles from (i-l)xlO° to i x 10°. We define the momentum 
flows, Xi{i = 1,9), into the i-th interval as a scalar sum of the momenta of all charged tracks and neutral showers 
directed in that interval. The momentum flows in corresponding forward and backward intervals arc combined 

Angular momentum conservation provides some additional discrimination between BB and continuum qq events. 
In qq production, the direction of the candidate thrust axis with respect to the beam axis in the cm. frame, 9t, tends 
to reproduce the 1 + cos^ 6t distribution of the primary quarks. The direction of the B candidate thrust axis for BB 
events is uniform. The B candidate direction with respect to the beam axis, 9b, exhibits a sin^ 6b distribution for 
BB events and is uniform for qq events. 

A Fisher discriminant is formed from 11 variables: the nine momentum flow variables, | cos 6*^1, and | cos 6*3!. 
The discriminant, J-, is the linear combination 

11 

1=1 

of the input variables, Xi, that maximizes the separation between signal and background. The coefhcients ai are 
determined from Monte Carlo simulation using a large set of continuum events and signal events modeled as B~^ — > 
if+TT+TT". We use the same set of coefficients ai for all three-body final states. The separation between the mean 
values of the signal and background distributions is approximately 1.3 times the signal width. 

For the A'tttt and KKt: final states, we impose a requirement on the Fisher discriminant variable that rejects 
90% of the remaining continuum background with about 54% efficiency for the signal. For the KKK final state, the 
continuum background is much smaller and we make a looser requirement that rejects 53% of continuum background 
with 89% efficiency for the signal. 

To determine the dominant sources of background from other decay modes of B mesons, we use a large set of 
Monte Carlo generated BB events where both B mesons decay generically Q. Most of the BB related background 
is found to originate from 5+ Z)°7r+, B+ J/ipK'^ and B~^ ^{2S)K^ decays. To suppress this type of 
background we apply the requirements on the invariant masses of the two-particle combinations that are described 
below. The background from B semileptonic decays is additionally suppressed by the electron veto requirement. The 
most significant background to the K^tt^t:'' final state from rare B decays is found to originate from B'^ rj' 
followed by ry' — ^ 7r+7r~7. We expect about 3% of these events to satisfy all the selection criteria. We find no significant 
background to the K'^K'^K^ final state from other known rare decays of B mesons. 

V. RESULTS OF THE ANALYSIS 

A. B+ K+TT+TT- 

For B^ AT+TT+TT^ decays, we form B candidates from three charged tracks where one track is positively identified 
as a kaon and the other two tracks are consistent with a pion hypothesis. Figure ^ shows the Dalitz plot for 
selected i?+ K^tt^tt^ candidates in the B signal region. Large contributions from the _B+ — > Z)°7r+, 
K+TT- and B+ J/ip{ip{2S))K+ , J/ip{ip{2S)) n+n' are apparent in the Dalitz plot. The J/-0(V'(25)) modes 
contribute to this final state due to muon-pion misidentification; the contribution from the J/'ip{'ip{2S)) — > e^e~ 
submode is found to be negligible (less than 0.5%) after the electron veto requirement. For further analysis, we 
exclude and J/ip{ip{2S)) signals by imposing requirements on the invariant masses of two intermediate particles: 
\MiK+n-) - Md\ > 0.10 GeV/c2; \M{h+h-) - Mj/^\ > 0.07 GeV/c^; \M{h+h-) - M^(2S)| > 0.05 GeV/c^, 
where /i+ and h~ are pion candidates. For the J /il){'ip{2S)) rejection, we use the muon mass hypothesis for charged 
tracks to calculate M{h^h^). To suppress the background caused by n/K misidentification, we exclude candidates 
if the invariant mass of any pair of oppositely charged tracks from the B candidate is consistent with the D Kit 
hypothesis within 15 MeV/c^ 2.5a), independently of the particle identification information {D veto). The Ai? 
and il/bc distributions for the remaining events are presented in Figs. §(a) and §(b), respectively. Here a significant 
enhancement in the B signal region is observed; the result of a fit to the AE distribution is presented in Table |. 
The expected AE and Afbc background distributions, which are the sum of luminosity-scaled below-threshold data 
and generic BB Monte Carlo events, are shown as open histograms in Figs, ^a) and|2|(b); the contributions from the 
BB Monte Carlo sample are shown as hatched histograms. There are no three-body charmless decays included in the 
generic BB Monte Carlo. 

To examine possible intermediate two-body states in the observed B'^ K~^'!T~^Tr~ signal, we analyze the K~^Tr~ 
and Tr~^TT~ invariant mass spectra shown in Fig. 0. To suppress the feed-across between the tt+tt" and if+Tr" states 
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we require the ii'+Tr" (tt+tt") invariant mass to be larger than 2.0 (1.5) GeV/c^ when making the tt+tt" {K~^Tr~) 
projection. The hatched histograms shown in Fig. ^ are the corresponding two-particle invariant mass spectra for the 
background events in the AE sidebands plotted with a weight of 0.62. 

The K~^Tr~ invariant mass spectrum is characterized by a narrow peak around 0.9 GeV/c^ which is identified as 
the K*{892)'^ and a broad enhancement above 1.0 GcV/c^ which is subsequently referred to as A'x(1400). In the 
TT+TT^ invariant mass spectrum two distinct structures in the low mass region are observed. One is slightly below 
1.0 GeV/c^ and is identified as the /o(980) while the other is between 1.0 GeV/c^ and 1.5 GeV/c^ and is referred to as 
/x(1300). Some excess of signal events can be also observed in the p°(770) mass region. The peak around 3.4 GeV/c^ 
is consistent with the process XcoK~^, Xco tt+tt^, and is the subject of a separate analysis Q. In this paper 

we exclude the —^ XcoK'^ candidates from the analysis of two-body final states by applying the requirement on 
the TT+TT" invariant mass: |M(7r+7r") - Af^^J > 0.05 GeV/c^. 

For further analysis we subdivide the full Dalitz plot area into seven non-overlapping regions as defined in Table ||. 
Regions from I to V are arranged to contain the major part of the signal from the B'^ i^*(892)°7r+, 5+ —* 
Js:x(1400)7r+, B+ p^{770)K+, B+ foi980)K+, and B+ fx{1300)K+ final states, respectively. The area in 
the Dalitz plot where Kn and tttt resonances overlap is covered by region VI, and region VII covers the rest of the 
Dalitz plot. The results of the fits to the AE distributions for all seven regions are summarized in Table ^ The 
procedure used for the extraction of the two-body branching fractions is described in detail in Section VI. 

B. B+ K+K+R- 

For the selection of 5+ K^K^K~ events, we use combinations of three charged tracks that are positively 
identified as kaons. The Dalitz plot for selected B^ — > K'^K'^K~ candidate events in the B signal region after the D 
veto is shown in Fig. ^. Since in this case there are two same-charge kaons, we distinguish the K~ combinations 
with smaller, M{K^ K')^i^, and larger, M(if +if ~)max, invariant masses. We avoid double entries per candidate 
by forming the Dalitz plot as Af ^(i4r+i4r~)max versus M^{K'^ K~)^in- The signal from the Cabibbo-suppressed 
— > D'^K^, K^_K~ decay mode is apparent as a vertical strip in the Dalitz plot. The corresponding 

Cabibbo-allowed i?+ L)"7r+, K^K^ decays can also contribute to this final state as a result of pion- 

kaon misidentification. We exclude candidates consistent with the i?+ D'^h^ hypothesis from further analysis by 
imposing the requirement on the K^K~ invariant mass \M{K^ K^) — Mjjo\ > 0.025 GeV/c^. The AE and Afbc 
distributions after the exclusion of D mesons are presented in Figs. ||(c) and ^(d) respectively. A large peak in the B 
signal region is apparent in both distributions. The result of a fit to the AE distribution is presented in Table 0. 

The K^K" invariant mass spectra for events from the B signal region are shown as open histograms in Figs.]5|(a)- 
H(c). The hatched histograms show the corresponding spectra for background events in the AE sidebands, plotted 
with a weight of 0.62. The M{K^K~),^in spectrum, shown in Fig. ||(a), is characterized by a narrow peak at 
1.02 GeV/c^ corresponding to the (?!)(1020) meson and a broad structure around 1.5 GeV/c^; this is subsequently 
referred to as /x(1500). To plot the M{K^K~)niax mass spectrum we subdivide the Af(iir+i4r~),„i„ mass region into 
two ranges: M{K+K~)^in < 1.1 GeV)c^ and M{K+K-)^in > 1.1 GeV/c^. The M{K+K-)^a^ mass spectra for 
these two regions are presented in Fig. [^(b) and Fig. ^(c) respectively. The prominent structure observed in Fig. ^(b) 
reflects the 100% (j) meson polarization in the B'^ (pK^ decay due to angular momentum conservation. In contrast, 
the distribution of signal events in Fig. |5|(c) is quite uniform after the background is subtracted. For the analysis of 
two-body final states we exclude events that are consistent with the i?+ XcoK~^, XcO K^K^ decay by applying 
the requirement on the K^K~ invariant mass: \M{K~^ K~) — Af^^ol > 0.05 GeV/c^. 



For further analysis we subdivide the full Dalitz plot area into the four non-overlapping regions defined in Tabic III . 
Regions I and II are arranged to contain the major part of the signal from the i?+ 0(lO2O)ii'+ and _B+ 
/x(1500)ii'+ final states, respectively. Regions III and IV cover the remaining part of the Dalitz plot. The results of 



the fits to the AE distributions for all four regions are summarized in Table III 



C. B+ -> K-TI+TI+, B+ K+K+TT- and B+ -> K+K-tv+ 

In general, we do not expect any signal in the B+ K^tt^tt^ and B^ K^K^tt^ final states. The Standard 
Model prediction for the i?+ K^K^tt^ branching fraction is of the order of 10^^^, and even much smaller for the 
5+ i^~7r+7r+ final state |l|]. However, these signals could be significantly enhanced in some extensions of the 
Standard Model |l2j, and, thus, these modes can be used to search for physics beyond tjie Standard Model. 

For the K'^K'^tt' final state, we reject candidates that are consistent with the i?+ D^K^, — ^ K^n~ decay by 
imposing the requirement on the K^tt~ invariant mass \M{K~^tt~) — Mdo\_> 0.10 GeV/c^. In case of the K^K^tt'^ 
channel we reject candidates that are consistent with the B~^ — > D°7r+, D'^ — K~^K~ decay with the requirement 
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\M{K^ K~) — M]jo\ > 0.05 GeV/c^. We also apply the D veto requirement for the three modes. The resulting 
AE and Mbc distributions for the K^tt^tt~^, K^K~^'k~ and K~^K~'k^ final states are presented in Figs. §(e)-||(j). 
Although we do not observe any signal in the AE distributions, there is an excess of events in the signal region of 
the Afbc distributions for the iir~7r+7r+ and K^K^tt'^ final states. These excesses could be caused by incorrectly 
reconstructed B decays. To subtract this background we subdivide the AE region into ten bins of 40 MeV width and 
determine the signal yield in each bin from the fit to the corresponding Mbc spectrum. The results of the fit, along 
with the expected contributions from the generic BB decays and the feed-down due to the particle misidentification 
from the B^ X+tt+tt" and B'^ K^K~ decay modes, are presented in Fig. ^. The latter two components 

are shown in Fig. || by the dotted and dashed histograms, respectively. The excess of events over the total expected 
background in the AE' signal region (two bins around AE = 0) is considered to be a signal yield. The results are 
summarized in Table ||. We do not observe a statistically significant signal in any of these three-body modes. 

The feed-across between K^tt^tt~ and K~^K~^K^ final states is found to be negligible. True B~^ K'^K~^K~ 
events reconstructed as K~^7t~^7t~ contribute mainly to the AE < — O.fO GeV region that is excluded from the fit. 
The fraction of true if+7r+7r~ events improperly reconstructed as K^K^K^ is less than 0.1%. 



VI. BRANCHING FRACTIONS 



To determine branching fractions, we normalize our results to the observed B^ — s- D^ir^, D'^ K^tt^ signal. 
This removes systematic effects in the particle identification efficiency, charged track reconstruction efficiency and the 
systematic uncertainty due to the cuts on event shape variables. We calculate the branching fraction for B meson 
decay to a particular final state / via the relation 

B{B+ ^ /) = B{B+ ^ D\+)B{D'^ ^ K+TT-)^ — , 

where Nf and Ndt^ are the numbers of reconstructed events for the final state / and that for the reference process, 
respectively; £/ and edtt are the corresponding reconstruction efficiencies. We use the signal yield extracted from the 
fit to the corresponding AE distribution; we do not use the Mbc distribution, because it, in general, suffers more from 
the BB background. 

The number of signal events for the reference process B^ D'^tt^ , — > if "•"tt" is found to be 1349 ± 40 for the 
K~^K~^K~ selection requirements and 805 ± 32 for the requirements used for all other three-body combinations. The 
corresponding reconstruction efficiencies are 26.8% and 16.1%, respectively. The reconstruction efficiency for each 
three-body final state is determined from the Monte Carlo simulation of events that are generated with a uniform 
Dalitzplot distribution. The branching fraction results for K^tt^tt" and K^K^K~ final states are presented in 
Table |, where the first quoted error is statistical and the second is systematic. The dominant sources of systematic 
error are listed in Table |^ We estimate the systematic uncertainty due to variations of reconstruction efRciency 
over the Dalitz plot using two sets of MC data generated with uniform distribution (phase space) and using some 
model (described below) . The uncertainty due to the particle identification is estimated using pure samples of kaons 
and pions from K~^tt~ decays, where the D° flavor is tagged using £)°7r+ decays. To estimate the 

uncertainty due to the signal and background shapes parameterization, we fit the AE distributions using different 
functions for the background description (linear, parabolic, exponential plus constant) and varying the parameters of 
the signal function (sum of two Gaussians with the same mean) within their errors. 

Since we do not observe a statistically significant signal in the K^tt^tt^, K^K^tt^ or K^K^tt^ final states, we 
place the 90% confidence level upper limits on their branching fractions. These limits are given in Table ||. To calculate 
the upper limits, we follow the PDG recommendation p^. 



A. Exclusive two-body branching fractions in the K tv n final state 

In the determination of the branching fractions for exclusive two-body final states, we have to take into account the 
possibility of interference between wide resonances. This requires some assumptions about the states that are being 
observed and, as a consequence, introduces some model dependence into the extraction of the exclusive branching 
fractions. The present level of statistics does not permit unambiguous interpretation of the Kx{i^OO) and /x(1300) 
states and, thus, it is not possible to use the data to fix all of the input model parameters. For this analysis we 
assume that the observed -firx(1400) and /x(1300) states are 0+ scalars. While this does not contradict the observed 
signal, some contributions from vector (1^) and tensor (2+) resonances cannot be excluded. The uncertainty related 
to this assumption is included in the model-dependent error described below. We ascribe to the Kx (14:00) state the 
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parameters of K^{U30) [M = 1412 MeV/c^, T = 294 MeV) and to /x(1300) state the parameters of /o(1370) (M = 
1370 MeV/c2, T = 400 MeV) 0. 

For further analysis we make the following assumptions: 

• The observed signal in the K^-k^tt^ final state can be described by some number of two-body final states. We 
restrict ourselves to the following set of exclusive two-body final states: K*(892)''7r+, iirx(1400)7r+, p°{n{))K+, 
/o(980)i^+ and /x(1300)iir+. We enumerate these final states as 1 through 5 in the order mentioned above. 

• Given this set of two-body final states, we determine the exclusive branching fractions neglecting the effects of 
interference. The uncertainty due to possible interference between different intermediate states is included in 
the final result as a model-dependent error. 

In order to extract the signal yield for each two-body final state, we perform a simultaneous likelihood fit to the 
/S.E distributions for the seven regions of the K^-k'^tt~ Dalitz plot. We express the expected number rij of signal 
events in the j-th region of the Dalitz plot as a linear combination 

5 

1=1 

where Ni is the total number of signal events in the j-th two-body final state and £ij is the probability for the i-th 
final state to contribute to the j-th region in the Dalitz plot. The Sij matrix is determined from the Monte Carlo 
simulation and includes the reconstruction efficiency. This procedure takes into account the effect of correlations 
between different channels in the determination of the statistical errors. 

The results of the fit are summarized in Table Combining all the relevant numbers, we calculate the product of 
branching fractions B{B'^ Rh^) x B{R h^h^), where R denotes the two-body intermediate resonant state. We 
present three types of errors for the branching fractions: the first error is statistical, the second is systematic, and the 
third refiects the model-dependent uncertainty. In general, the model-dependent error is due to uncertainties in the 
effects of interference between different resonant states. We estimate this error by means of a _B+ if +7r+7r~ Monte 
Carlo simulation that includes interference effects between all the final states mentioned above. We vary the relative 
phases of the resonances and determine the signal yield using the procedure described above. The maximal deviations 
from the central values are used as an estimate of the model dependence of the obtained branching fractions. We 
find that the model-dependent errors associated with the wide resonances introduce significant uncertainties into the 
branching fraction determination. In the case of the p'~'{770)K^ final state, this effect is enhanced by the smallness 
of the signal itself. Since we do not observe a significant signal in this channel, we report a 90% confidence level 
u pper limit. The s tatistical significance of the signal, in terms of the number of standard deviations, is calculated as 
■y/— 2 ln(>Co /'Cmax ) 7 where /!niax and Cq denote the maximum likelihood with the nominal signal yield and with the 
signal yield fixed at zero, respectively. 

B. Exclusive two-body branching fractions in the K'^ final state 

In the case of the three charged kaon final state, we clearly observe the (/)(1020) meson plus a very broad /x(1500) 
structure that we currently cannot interpret unambiguously. It could be a complex superposition of several interme- 
diate states and some contribution from the non-resonant i?+ K^K'^K' decay is also possible. For our study of 
systematic and model-dependent uncertainties, we construct a simplified model and parameterize the /x(1500) struc- 
ture as a hypothetical scalar state with M — 1500 MeV/c^ and F = 700 MeV. We find qualitative agreement between 
the experimental Dalitz plot distribution of the signal events and that obtained from the Monte Carlo simulation with 
this simple model. 

Then we extract the signal yield for the two-body final states; _B+ (f){lQ20)K^ and the so-called 5+ 
fx{f500)K~^, which is, in fact, all of the remaining signal. We follow the same procedure as we used for the K^tt^tt^ 
final state. The signal yields are determined from a simultaneous fit to the AE distributions for four separate regions 
of the K^K^ K~ Dalitz plot. The results of the fit are summarized in Table [V^ . 

We determine the model-dependent error in the same way as we did for the X+tt+tt^ final state. In the case of the 
K'^K^K~ final state the model-dependent error is found to be much smaller than in the K^tt^tt^ final state. This 
is mainly due to the small width of the 0(1020) meson. 
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VII. DISCUSSION & CONCLUSION 



The high quaUty of tt/K separation at Belle allows us to measure, for the first time, the branching ratios for the 
three-body modes B{B+ -f K+n-n+) = (55.6±5.8±7.7) x 10"^ and B{B+ K+K-K+) = (35.3±3.7±4.5) x IQ-^ 
without assumptions about particular intermediate mechanisms. CLEO and BaBar p5| have previously placed 
upper limits on the branching fractions of non-resonant three-body decays: B{B^ K^n'^^) < 28 x 10~^ (CLEO), 
B{B+ K+n+TT-) < 66 X IQ-*^ (BaBar), B{B+ K+K+R-) < 38 x IQ-^ (CLEO). A comparison of the applied 
selection criteria shows that CLEO and BaBar restricted their analyses to the region of invariant mass above 2 GeY /(? 
for any pair of the particles. This requirement effectively removes most of the low mass resonances that provide the 
dominant contribution to our observed signal. They assume a uniform distribution of events over the Dalitz plot to 
obtain the limits quoted above. The upper limits reported here for the if^7r+7r+, K^K^tt" and K'^K~'k'^ modes 
are considerably more restrictive than previous limits from CLEO [Q and OPAL . 

Significant signals are observed for the first time in the decay modes B+ /o(980)i4:+ and B+ K*{892fTT+. The 
measured branching fraction product for the fo{980)K+ final state is B{B+ fo{980)K+) x ;B(/o(980) tt+tt") = 
(9.6t2'3li'5lo g) ^ 10^^. This is the first observation of a -B decay to a charmless scalar-pseudoscalar final state. 
The measured branching fraction product for the i^*(892)°7r+ final state is B{B+ K*{892fTT+) x B{K*{892f 
K+-K-) = il2.9tlil\A-l^5) X 10-6. Using the value of B{K*{892)° K+n-) = 2/3, we translate our measurement 
into the branching fraction B{B+ X*(892)°7r+) = [19 Atfltl'^tll) x 10"^. The significant enhancement in the 
K~^TT~ invariant mass spectrum above the K*{892) mass agrees with the scalar ^0(1430) hypothesis. This is also in 
agreement with theoretical predictions for the B'^ iiTg (1430)7r+ branching fraction based on the factorization 
model. Nevertheless, we cannot exclude some contribution from the tensor i4r|(1430) state. 

The interpretation of the peak with a tt+tt" invariant mass around 1300 MeV/c^ in the X+tt+tt^ system is less 
certain. There are two known candidate states: the /2(1270) and /o(1370) ||l^. Attributing the peak to the /o(1370), 
with its rather small coupling to tt~^'!t~ JlSt , would lead to an unusually large branching fraction for a charmless B 
decay mode. On the other hand, as recently shown in | p9[ , the factorization model predicts a very small branching 
fraction for the /2(1270)iir+. If our observation is, in fact, due to the /2(1270), this would provide evidence 

for a significant nonfactorizable contribution. 

We cannot identify the broad structure observed in the 5+ K^K^K^ final state above the (/)(1020) meson. It is 
hardly compatible with the presence of a single scalar state, either /o(1370) or /o(1500) @. We also cannot exclude 
the presence of a non-resonant contribution or the case of several resonances contributing to the excess in the K'^ K~ 
invariant mass spectrum seen around 1.5 GeV/c^. 

We find that effects of interference between different two-body intermediate states can have significant influence on 
the observed two-particle mass spectra and a full amplitude analysis of three-body B meson decays is required for a 
more complete understanding. This will be possible with increased statistics. 
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TABLE I: Branching fractions and 90% C.L. upper limits for S"*" — » K'^K^h final states. 



Three-body Efficiency 


Yield 


B 


mode 


(%) 


(events) 


(io-«) 




17.3 


237 ± 23 


55.6 ±5.8 ±7.7 


K+K+K- 


24.0 


210 ±21 


35.3 ±3.7 ±4.5 


K--K+H+ 


16.2 


12 ±9 


< 7.0 


K+K+n- 


14.2 


2.0 ± 5.3 


< 3.2 


K+K-n+ 


14.6 


26 ±12 


< 12 



TABLE II: Results of the fit to the AE distribution for different regions in the K^-k^tt Dalitz plot. Columns list the definition 
of each region, reconstruction efficiency from Monte Carlo simulation and signal yield. 



Dalitz plot 


Mass range 


Efficiency Yield 


region 


(GeV/c^) 


(%) (events) 


I 


M{Ktv) < 1.00; M(7r7r) > 1.50 


20.7 ±3.7 47.11^:? 


II 


1.00 < M{KTr) < 2.00; M{Tnv) > 1.50 


19.2 ± 1.3 56.0t\l 


III 


M{-ktt) < 0.90; M{KTr) > 2.00 


16.7 ±2.4 17.7tr:t 


IV 


0.90 < MiwTT) < 1.06; M{Ktt) > 2.00 


19.9 ± 3.2 34.7I7 3 


V 


1.06 < M(7r7r) < 1.50; M{K7r) > 2.00 19.6 ± 1.7 33.4l|^ 


VI 


M{K7v) < 2.00; M{tt7v) < 1.50 


14.7 ±3.3 14.9l^;? 


VII 


M{Kn) > 2.00; M{mr) > 1.50 


16.1 ±0.5 ll.Slgj 



TABLE III: Results of the fit to the AE distribution for different regions in the K^K'^K Dalitz plot. Columns list the 
definition of each region, reconstruction efficiency from Monte Carlo simulation and signal yield. 



Dalitz plot 
region 


Mass range 
(GeV/c^) 


Efficiency 
(%) 


Yield 
(events) 


I 

II 


M{KK)min < 1.04 

1.04 < M{KK)n,in < 2.00 


24.6 ± 2.5 
23.3 ± 0.8 





III M{KK)min > 2.00; M{KK)^^^ > 3.40 23.9 ± 1.1 

IV ..l/(A'iv)nnu > 2.00: M(KK < 2 1.7 ±0.8 


:i2.:r;::« 




TABLE IV: List of systematic errors (in percent) for the B'^ 


K+h+h- 


branching 


fractions. 


Source K'^-k'^it~ K~7t'^-r 


+ K+K-T,+ 


K+K+t:- 


K+K+K- 


B —f Dtt and D Kit branching fractions 9.7 9.7 
efficiency nonuniformity over the Dalitz plot 7.6 
background and signal parameterization 6.3 
particle identification 


9.7 
3.0 


9.7 
3.0 


9.7 
3.7 
4.7 
6.0 


total 13.8 9.7 


10.1 


10.1 


12.9 
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TABLE V: Results of the simultaneous fit to the K'^tt'^it final state. 



Two-body 


Efficiency 


Yield 


Significance B 


B+^Rh+ ^ ^R^h+h- 


mode 


(%) 


(events) 


(^) 


(io-«) 


A'* (892)%+ 


18.9 




6.2 


1 9 Q+2.8+1.4+2.3 
-■-^•^-2.6-1.4-4.5 


Kx{UOO)tt+ 


16.2 




4.9 


-lA r+3.5+1.8+3.3 
-■■*-^-3.3-1.8-6.B 


p°{770)K+ 


15.1 




0.8 


< 12 


M980)K+ 


17.8 




5.0 


Q „+2.5+1.5+3.4 
^•"-2.3-1.5-0.8 


fx{1300)K+ 


16.9 


A6t\t 


3.9 


11 1+3.4+1.4+7.2 
-'^-'^•^-3.1-1.4-2.9 



TABLE VI: Results of the simultaneous fit to the K+K+K' final state. 



Two-body 


Efficiency 


Yield 


Significance Bb+^rh^ x B^^^+h.- 


mode 


(%) 


(events) 


(a) (10-«) 


(t>{W20)K+ 


23.6 


421?;^ 


Y 9 7 9+1.5+0.9+0.4 
'•^ '•^-1.4-0.9-0.4 


fx{1500)K+ 


21.3 


14611^ 


19 97-^+3.2+3.5+1.4 
-'^ ^'•"-3.2-3.5-1.4 
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FIG. 1: The Dalitz plot for — > it'+TT+Tr candidates from the B signal region. 
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FIG. 2: The AE (left) and Mbc (right) distributions for B+ K+h+h~ final states: (a, b) - AT+tt+tt"; (c, d) - K+K+K'; 
(e, f) - K~-w'^n'^; (g, h) - K'^K'^tt~; (i, j) - K'^ K~ ir'^ . Points with errors represent data, open fiistograms arc the proper sum 
of the bclow-threshold data and BB Monte Carlo; the hatched histograms show the contribution of BB Monte Carlo only. The 
solid lines display the signal plus background combined shape. The dashed lines correspond to the background shape only. 
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FIG. 3: The (a) K^-k and (b) ■k'^-k invariant mass spectra for selected -B"*" — > K'^n'^n candidates in the B signal region 
(open histograms) and for background events in the AE sidebands (hatched histograms). 
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FIG. 4: The Dalitz plot for — > JsT candidates from the B signal region. 
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FIG. 5: The K'^ K~ invariant mass spectra for selected — > K'^ K~ candidates in the B signal region (open histograms) 
and for background events in the sidebands (hatched histograms), (a) The K'^K~ combination with the smaller invariant 
mass. The inset shows the 0(1020) region in 2 MeV/c^ bins; (b) the M{K+K-)^s^ spectrum with M{K+K-)^i^ < 1.1 GeV/c^ 
and (c) the M{K+K-)^^^ spectrum with M{K+K-)^in > 1-1 GeV/c^. 
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FIG. 6: Results of the fit to the Mbc distributions in AE bins. The data (points with errors) are compared with the MC 
expectation (open histogram). The feed-down from K'^K'^K~ and K'^'K'^n~ final states is shown by the dotted and dashed 
histograms, respectively. 



